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o : 

f*^ , The statistical mechanics of a ribbon polymer made up of two semiflexible chains is studied 

■ using both analytical techniques and simulation. The system is found to have a crossover transition 

' at some finite temperature, from a type of short range order to a fundamentally different sort of 

short range order. In the high temperature regime, the 2-point correlation functions of the object 
are identical to worm-like chains, while in the low temperature regime they are different due to a 
twist structure. The crossover happens when the persistence length of individual strands becomes 
comparable to the thickness of the ribbon. In the low temperature regime, the ribbon is observed to 
have a novel "kink-rod" structure with a mutual exclusion of twist and bend in contrast to smooth 
worm-like chain behaviour. This is due to its anisotropic rigidity and corresponds to an infinitely 
strong twist-bend coupling. The double-stranded polymer is also studied in a confined geometry. It 
is shown that when the polymer is restricted in a particular direction to a size less than the bare 
^ . persistence length of the individual strands, it develops zigzag conformations which are indicated 

by an oscillatory tangent-tangent correlation function in the direction of confinement. Increasing 
the separation of the confining plates leads to a crossover to the free behaviour, which takes place 
at separations close to the bare persistence length. These results are expected to be relevant for 
experiments which involve complexation of two or more stiff or semiflexible polymers. 
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I. INTRODUCTION AND SUMMARY 

> 

' There has been a lot of recent interest in physical properties of biopolymers, ranging from elasticity of biopolymer 
networks and its use in the prediction of mechanical properties of cells to direct visualisation of single chain properties. 
Examples of important biological macromolcculcs whose physical properties have been recently studied are: actin, a 
double-stranded semiflexible protein polymer which forms an integral part of the cytoskeleton (mechanical structure) 
of eukaryotic (e.g. fungi, plants, animals) cells; microtubules, multi-stranded rigid and dynamic protein polymers 
which form one of the main components of the cytoskeleton of eukaryotic cells and play an important part in their 
organisation; DNA, which carries the genetic code of all living organisms 

Since many of the processes involved in cell function(e.g. DNA replication in cell division) require major structural 
""j" 1 , changes of these biopolymers , there is a need for more microscopic but still analytically tractable models of such 
polymers which go beyond the simple picture of such molecules as homogenous elastic rods ||. Motivated by this, we 
study such a microscopic model specifically looking for qualitative differences between the behaviour of such molecules 
and simple worm-like chains [Q. In addition, most analyses of the worm-like chain models of polymers have focused 
on ground state properties (long chains) or bulk quantities ||. It is interesting to look at the effect of fluctuations, 
1 spatial correlations and finite size on these systems. 

A double-stranded semi-flexible polymer chain is the basic structure of many biopolymers. Examples are of double- 
stranded biopolymers are DNA and proteins such as actin. The model most used in the study of biopolymers is 
that of the worm-like chain Q in which the polymer flexibility (structure) is determined by a single length, the 
persistence length £ p which measures the tangent-tangent correlations. For example, DNA has a persistence length 
lp « 50nm whilst for actin £ p « 17/im. These biopolymers are known to have a more complex 'twisted' structure. 
The multi-stranded nature of these polymers is also not taken into account in a simple worm-like chain model. It 
is not clear if such a fine structure will have an effect on the global properties of these objects. A possible effect of 
such fine structure is what we attempt to study in this article. Our model is, in a sense, microscopic because the 
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interaction between the bend and twist degrees of freedom is a result. This is fundamentally different from previous 
approaches Jjl which try to include the twist degrees of freedom by adding extra terms to the free energy. 

In a previous letter Q we studied a version of the rail-way track model of Everaers-Bundschuh-Kremer (EBK) || 
for a double-stranded semi- flexible polymer, embedded in a d-dimensional space for arbitrary d. The main purpose 
of this article is to present a detailed description of our theoretical and numerical calculations and in addition we 
present some new results on the effect of confinement on the statistical mechanics of the ribbon polymer || . Excluded 
volume and electrostatic interactions have been ignored throughout. 

We find that the system, has qualitatively different properties in the low temperature and high temperature regimes, 
in contrast to what one might naively expect from an inherently one-dimensional system with local interactions and 
constraints. The tangent-tangent correlation function decays exponentially in the whole range of temperatures with 
a "tangent-persistence length" ^tp that has a very slow temperature dependence, and whose scale is determined by 
the (bare) persistence length of a single strand £ p = n/ksT (k is the bending stiffness of a single strand). Note that 
it is independent of a, the separation of the two strands, which is the other relevant length scale in the problem. 
However, the correlation function of the bond-director field, defined as a vector that determines the separation and 
coupling of the two strands of the combined polymer system, has different behaviour below and above the temperature 
T x ~ 4.27n/dkBa. While it decays purely exponentially for T > T x , there are additional oscillatory modulations 
for T < T x . The related "bond-persistence length" £ BP does not change appreciably at high temperatures, where 
its scale is again set by £ p alone. In the low temperature phase, however, £bp does show a temperature dependence. 
In particular, £ BP ~ ^ /3 a 2 / 3 oc T" 1 / 3 for T ~ 0, while £ BP ~ l v for T ~ T x . Similarly, the "pitch" H, defined 
as the period of oscillations in the low temperature regime, changes drastically with temperature , ranging from 
H r~j £ BP ~ £p^ 3 a 2 / 3 near T = 0, to H ~ near T — T x . At T = we restore a flat ribbon which has true long-range 
order in both the tangent and bond-director fields. The ribbon is essentially a rigid rod. As we approach T = 0, the 

1/3 

persistence lengths and the pitch diverge with the scaling H ~ £ BP ~ £^ p . 

The spontaneous appearance of a short range twist structure may be understood in the language of the homogenous 
rod models Q as a local twist-bend coupling, which is observed up to a screening length £ BP that lies between a and 
£ p . We find that the anisotropy in the rigidity of a ribbon results in a "kink-rod" structure, in which the ribbon is 
at every point along its contour either twisted and unbent (rod) or bent and untwisted (kink). This inhomogeneous 
behaviour is in sharp contrast with uniform worm-like chain behaviour, and can be interpreted as an infinitely strong 
twist-bend coupling. This structure is, however, screened on long length scales due to the fact that the ribbon is a 
one dimensional system with short range interactions. The short range twist order and the kink-rod structure will 
naturally disappear when a ~ £ p corresponding to the temperature T x . We also observe a twist-stretch coupling. 

We also study the effects of confinement on the double-stranded polymer. We enforce the confinement to a box of 
size R as an additional constraint. For R 3> £ p , we recover the behaviour of free double-stranded semiflexible polymers 
0. In particular, the tangent-tangent correlation has a purely exponential decay with a characteristic length scale of 
order £ p , while the bond-director field develops a crossover to a phase with oscillatory correlations. As R is decreased, 
there is a crossover to a phase with oscillatory tangent-tangent correlations about R ~ £ p . For we find that 

both the persistence length in the perpendicular direction £±, and the characteristic oscillation length A, scale as 
(i'pi? 2 ) 1 / 3 which is considerably less than £ p . The same crossover in the bond-director field also persists in this limit, 
and in particular, there is a regime in which both the tangent and the bond-director correlations are oscillatory. 

We introduce and define our model in Sec. ||, and then describe a mean field approach in Sec. Ill, which can be 
used to obtain closed form expressions for various correlation functions. In Sec. |^ we discuss a physical argument 
for the results we obtain using a plaquette model. In Sec. ^ we describe some extensive Molecular Dynamics/Monte 
Carlo simulations, which we use to calculate the correlation functions, and compare them with the mean field results 
of Sec. [II. We discuss the thermal kink- rod structure of stiff ribbons in Sec. VI. and the effects of confinement on 



the conformations of ribbon polymers in Sec. VII. Finally, Sec. VIII summarises our results and the limitations of 
our approach. 



II. RAILWAY TRACK MODEL 



To study the effect of a double-structure on semiflexible polymers, we consider a version of the railway track model 
of Everaers-Bundschuh-Kremer (EBK) In our approach, we are able to consider polymers embedded in a d- 
dimensional space for arbitrary d. The system is composed of two semiflexible chains, each with rigidity n, whose 
embeddings in d-dimensional space are defined by i"i(s) and r2(s). The Hamiltonian of the system can be written as 
the sum of the Hamiltonians of two wormlike chains 
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n = - i as 



d 2 ri(a) 
ds 2 



d 2 r 2 (s) 
ds 2 



(1) 



We assume that the individual strands (that make up the double-stranded polymer) are inextensible: (dri/ds) 2 = 
(dr 2 /ds) 2 = 1. The ribbon structure is then enforced by having r 2 (s') separated from r^s) by a distance a, i.e. 
r 2(s') = r i( s ) + ab(s), where \s — s'\ can be non-zero but is small. We have defined a bond-director field b(s), which 
is a unit vector perpendicular to both strands (see Fig. |l|). The chains are assumed to have permanent bonds (such 
as hydrogen bonds) that are strong enough to keep the distance between the two strands constant. In Ref. Q, it 
is argued that the relevant constraint on the system would then require that the arclength mismatch between the 
two strands in a bent configuration should be very small. We can calculate the arclength mismatch for the bent 
configuration as As = |r 2 (s) — ri(s) + a b(s)|, where a is the separation of the strands. We impose the constraint as 
a hard one, namely, we set As = as opposed to Ref. ||. Physically this means we do not allow bends in the plane 
of the ribbon. These bends are not important in d > 2 because as we shall see the lower length scale will be set by 
the 'pitch' which will make the in-plane fluctuations of the ribbon irrelevant [g] . 

We can argue that the simplifying assumption does not change the behaviour of the system. If we impose a soft 
constraint as in Ref. || using an energy term like (k/2) J (As) 2 , we can see that the length I = (n/ka 2 ) 1 / 2 determines 
two different regimes; the interesting one being L I [L is the length of the chains). Hence, our hard constraint in 
fact only restricts us to the case of interest. 

We implement the constraint As = by introducing the "mid-curve" r(s): 



r 2 (s) 



r(s) + - b, 
r(s) - - b. 



(2) 



In terms of the tangent to the mid-curve t = dr/ds, which we call the tangent- director field, and the bond-director b, 
the Hamiltonian of the system can now be written as 



H 



subject to the exact (local) constraints 



ds 



dt(s) 
ds 



("if)"- 1 . 



a 

T 



d 2 b(s) 



b 2 = l, 



(3) 



(4) 



a db N 
(t±^)-b = 0, 

For a weakly bent ribbon, the Hamiltonian in Eq.(||) can be conveniently thought of as having two major contributions: 
a bending energy Tij, (the first term), and a twisting energy TL t (the second term). 



III. MEAN FIELD THEORY: FREE CHAINS 



It is well known that the statistical mechanics of semiflexible chains are difficult due to the constraint of inextensi- 
bility. Various approximation methods have been devised to tackle the problem. A successful scheme, that somehow 
manages to capture the crucial features of the problem, is to impose global (average) constraints rather than local 
(exact) ones [p|-^3|. This approximation is known to be good for calculating the average end-to-end length. It can be 
used for the probability distribution of the end-to-end length only if the persistence length is much less than the chain 
length so that the chain conformation can be considered isotropic. We can get good insight to the approximation 
scheme, by considering the fact that it corresponds to a saddle-point evaluation of the integrals over the Lagrange 
multipliers, which are introduced to implement the constraints Jl2]JT^ |. In this sense, it is known to be a "mean-field" 
approximation in spirit. One can then go further by considering the effect of fluctuations on this mean-field result. The 
power of this approach is that one can easily calculate quantities which turn out to be very difficult if the constraints 
are required to hold exactly || . 

To study the effect of fluctuations, we have performed a systematic 1/d-expansion (sketched in Appendix |A|) [ pT[ . 
We see that no divergent behaviour appears when we calculate the diagrams of the 2-point correlation functions. This 
means that the mean-field behaviour of these functions, at least, won't change due to fluctuations, although it does 
not preclude differences in higher order correlation functions. 
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With the above discussion as justification, we apply the same approximation scheme to our problem defined in 
Sec. |l| The local constraints in Eq.([|) are relaxed to global ones. This can be done by adding the corresponding 
"mass terms" to the Hamiltonian 



7i„ 



knT 



ds 



a dbV 

2 dJy 



a dbY 

2 dJy 



ca 



(5) 



db 



2 ds 



+ t - - — -b + — t + - — • b 



a db 



2 ds 



where 6, c, and e are dimensionless constants. The partition function is then given by 

H + H n 



Z[3,K] = J Vt(s) Pb(s)exp 



J ds [J(s) • t(s) + K(s) • b(* 



(6) 



We next determine the constants self consistently by demanding the constraints of Eq.(Q) to hold on average, where 
the thermal average is calculated by using the total Hamiltonian H + Ti m . Note that in choosing the above form, we 
have implemented the "label symmetry" of the chains, namely, that there is no difference between two chains. It is 
convenient to take the limit of an infinitely long chain and perform the functional integrals in momentum space. The 
Fourier transforms are defined A(q) = J ds exp{zgs}A(s). We have 



Z[J,K] = J Vt(q) 2?b( 9 )exp j-i J ^(t(-g), b(-«)) • M(q) 
The Gaussian integration can then be easily performed. It yields 

Z[3,K] = exp 

where 



t(g) 
b(g) 



dq 
2tt 



J t + K b 



' r ^ft-q),K(-q)).M-Hq) 



J(«) 

K(g) 



M(q) = 



2£ p q 2 +4b/£ p 2e/£ p 

2e/£ p e p a 2 q 4 /2 + ba 2 q 2 /e p + ca 2 /e p 
The averages are easily obtained from Z, for example 

S 2 InZ 



(7) 



(8) 



(9) 



(U(q) ■ bj{q')) = 



The next step is to demand self consistently that 



t(s)± 



a db(s) 
2 ds 



(b(s) 2 ) =1, 

2 W, 



(10) 



t(s) ± a - Ah[s) 



b(s)) =0, 



2 ds 

The self-consistency leads to the following set of equations for the constants b, c and e: 

1 , a 2 Vc. _ 1 
4V26 ^ 4dl|" d 



e = 



(11) 



The above equations, which are nonlinear and difficult to solve exactly, determine the behaviour of b and c as a 
function of u = a/£ p . We have solved them numerically in d — 3 and the solutions are given in Fig. |[ One can solve 
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Eg. (pi]) analytically in two limiting cases. For u <C 1 we find b = d 2 /32 , and c = (d/ y/2) 4 / 3 u~ s / 3 , whereas for u ^> 1 
we find b = d 2 /8, and c = 4/u 4 . In Fig. [|, the behaviour of b and c is plotted as a function of u. Note that it is 
proportional to T and can be viewed as a measure of temperature. 

We can then calculate the correlation functions. For the tangent -tangent correlation one obtains 



<t(«).t(0)> = 



4v/26 



exp 



2b 



whereas for the bond-director field one obtains 

exp (—(b — VW 



(b( S ).b(0)> = 



d£ 2 p 



2a 2 y/b 2 



a/2 



exp 



-{b + VE^c) 1 / 2 



(b-VB^c) 1 / 2 



(b + Vb 2 ^) 1 / 2 



(12) 



(13) 



The tangent -tangent correlation (Eq.(|12|)) is exactly what we obtain for a single worm-like chain, and implies uniform 



behaviour for all temperatures. 



3) on the other hand, indicates a change of behaviour at b — c for the 



bond-director correlation. The correlation is over-damped for b 2 > c (high temperatures), while it is under-damped 
(oscillatory) for b 2 < c (low temperatures). The interesting point b 2 = c happens for u c = 16/(1 + \/2) 3 / 2 d ~ 4.27/d, 

that leads to the value for T x quoted above (see Fig. ||). We also find a divergence in the specific heat, C'v = fj§ 
where F — —kgT log Z at T x . It should be noted that it is not a thermodynamic phase transition in the sense of 
long-range ordering and broken symmetry. It is a crossover that appears due to competing effects, and the transition 
is from a state with some short-range order to a state with a different short-range order. Similar phenomena have been 
observed in Ising-like spin systems with competing interactions p5[ and the crossover (transition) point corresponds 
to a type of 'Lifshitz point' for a 1-d system. 

It is interesting to study the bond-director correlation in the limiting case b 2 <C c, that corresponds to relatively 
low temperatures. Using the asymptotic forms for b and c, one obtains 



(b(*).b(0))=V2 



exp 



d 



A£ p a 2 



1/3 



1/3 

U p a 2 I 8 + 4 



(14) 



for very low temperatures. From the above expressions for the correlation functions, one can read off the persistence 
lengths £ TP ~ £ P and £ BP ~ (^a 2 ) 1 / 3 , and the pitch H ~ (^a 2 ) 1 / 3 . 

^From the tangent-tangent correlation function, we can calculate the end-to-end distance. It yields 



((r( S )-r(0)) 2 ) = f f d Sl d S2 <t(si)-t(s 2 )) 



(15) 



o Jo 

s ■ 



46 



1 



-V2bs 



which is similar to worm-like chains. It interpolates between the limiting behaviours of random walks (~ (d£ p /4b)s) 
for s ^> £ p to rods (~ (d/4^/2b)s 2 ) for s <C l v - However, it is interesting to note that there is a shrinking in the length 
of the rod by a factor of (d/4-v/26) 1//2 , which varies smoothly from 1 at a <C t p to l/v2 at a 3> £ p . This implies that 
a polymer made up of two inextensible strands is always "slightly extensible" at any finite temperature, due to the 
presence of twist fluctuations. This is exactly the twist- stretch coupling studied by various authors using homogenous 
elastic rod models || . Note that in our model this coupling is a result, as opposed to the elastic rod models in which 
it must be added by hand. A microscopic model proposed by O'Hern et al M|, which describes DNA as a stack of 
plates, also predicts a twist-stretch coupling. 

A simple scaling argument can account for ^tp and ^bp in the low temperature regime. Consider applying a uniform 
bend of radius of curvature A to a section of ribbon of length A without twisting it. The corresponding bending energy 
(calculated using Hb) is given by Et, ~ ksT£ p /\. We can then estimate £tp by finding the length A for which the 
bending energy Eh becomes comparable to fcsT. Similarly, applying a uniform twist per length 27r/A to a section 
of the ribbon of length A without bending will cost a twist energy E t = ksTipa 2 / X 3 (calculated using Tit)- The 
wave-length A at which the twist energy E t becomes comparable to fcsT gives £bp- 



IV. PLAQUETTE MODEL AND COMPETITION 



The nature of competition in our double-stranded polymer system can be understood using a plaquette model. We 
can coarse-grain the ribbon to a length-scale (£) where we can consider it to be made up of plaquettes which are 
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joined up to form a ribbon (see Fig. ||). We can then define effective coarse-grained bond B^ and tangent Tj director 
fields for each plaquette. 

The energy expression corresponding to bends comes from the product of the tangent-directors of the neighbouring 
plaquettes: 



(3H b 



I 



-i+l- 



(16) 



In a spin analogy, this corresponds to a classical Heisenberg ferromagnet in 1-d, and has no competition . 

If we choose I -C l v in the coarse- graining process, we can safely assume that the ribbon is rod- like (we freeze out 
the bending modes) and that it only has twist fluctuations. We may then write b(s) = ei cos 6>(s) + 62 smd(s) where 
the (fixed) unit vectors ej, {i = 1,2} span the plane perpendicular to the rod. Rewriting Eq.(||) in terms of 6(s) and 
implementing the constraints, we obtain 



(3H 



La 2 



ds 



a 2 \d 2 0\ 2 \d s 6\< 
4-a 2 \d s 9\ 2 



(17) 



subject to the constraint: |d s 0(s)| < 2/ a. Note that the lowest order contribution to the twist potential starts from a 
quartic term. 

We can now expand the nonlinear term in the above Hamiltonian, and perform the coarse-graining, in the framework 
of a perturbation theory, by integrating out the modes between 1 /£ and \ j a in the momentum shell. We can then 
determine the form of the coarse-grained Hamiltonian and calculate the renormalizcd coupling constants. We keep 
only terms up to second (Gaussian) order, which is a good approximation for I < (ipa 2 ) 1 ^ 3 , and obtain 



PHt 



B 



i+l 



J 2 ^B, 



B 



i+2 



(18) 



where 



J\ = Co + C\ 



ipO 2 



A J l P a 

and J-2 = —Ci — 



P 



where Cj, {i = 0, • • • , 2} are constants of order unity. In contrast to the bending energy, the effective twist energy is 
frustrated due to the opposite sign of J± and J2 |15|. In the spin analogy, this would correspond to a model with 
next nearest neighbour competing interactions similar to the so-called ANNNI model, which develops oscillations for 
certain values of the ratio — J2/J1 (of order 1) |L5|. This corresponds to I ~ (ipO 2 ) 1 ^ 3 , and we can thus account for 
the pitch H ~ (V 2 ) 1/3 - 

The competition is present only at nonzero temperatures, and is merely due to topological constraints of the ribbon. 
Another, more physical way of understanding this competition is to consider the interaction between two neighbouring 
twisted regions. It is easy to see that twists of opposite sign meeting at an edge tend to unwind (annihilate) each 
other |17|, while twists of the same sign are trapped when they meet; they do not annihilate each other and add up 
(see FigTTl). 



V. SIMULATION 



An intriguing feature of the behaviour of this model is that, although the ground state (T = 0) configuration of the 
system is a flat ribbon, and supports no twists, upon raising the temperature, a twisted structure with short range 
twist order develops. We have confirmed this by performing extensive Molecular Dynamics (MD) /Monte Carlo (MC) 
simulations of double-stranded scmiflcxiblc polymers. A bead-spring model with bending and stretching energies was 
used. We combined a velocity Verlet MD coupled to a heat bath with an off-lattice pivot MC algorithm. The MD 
was useful for equilibrating the shorter length-scales and MC for the long length-scales. 

We used a triangular lattice to discretise the ribbon (see Fig. [|). The position of the ith bead is r 4 and we assume 
all the beads have mass m. The two chains making up the double strands join the odd ({1, 3, 5, ... , 799}) and even 
({2, 4, 6, ... , 800})beads together. The Potential Energy is given by 

^J^l = V fc,[(r <+ i - r,) 2 - i 2 ] + k s [(r i+2 - r,) 2 - *§] - h cos^, (19) 



G 



where 



cos 6^ = (r i+2 - ri) ■ (r; - ri_ 2 )- 

We have a bending constant kf, only for the springs joining beads on the same chain and a stretching constant k s for 
every spring. We also have a short-range repulsion between nearest neighbour beads. The MD simulation is performed 
by integrating a Langevin equation for every bead 

m^+r^ = -V n [/ + /,, (20) 

where fi is a random number chosen from a range set by T representing the heat bath. The simulations were performed 
at ksT = 1. We are in the dissipative regime so we can ignore the inertial term. The friction term is set to (T = 0.7) 
and the noise is chosen so as to satisfy the fluctuation dissipation theorem. The equilibrium bond length was set to 
£o = 1.6. The simulations were done with k s — 1000. 

We performed in general 10 6 integration time steps followed by 10 attempted pivot moves. A pivot move is 
an attempt to rotate a portion of the chain by a small random angle around a randomly chosen bead. The MC 

part is done with the usual metropolis algorithm accepting pivot moves with a probability exp ( ^ ■ This mixed 

MD/MC procedure was repeated 10 3 times until the configurations were equilibrated. Equilibration was checked by 
starting from crumpled chains and fully extended chains and verifying that the same values for radius of gyration 
and correlation functions was obtained. We simulated double-stranded ribbon chains of 2 x 400 monomers. The 
simulations were performed on a CRAY T3D with 128 processors allowing us to simulate 128 chains in parallel. 

Typical equilibrated polymer configurations, shown here in Figs. ^, |^ 0, suggest that at low temperatures the 
polymer can be viewed as a collection of long, twisted (straight) rods that are connected by short, highly curved 
sections of chain which we call "kinks" , as opposed to a smooth worm-like conformation. This structure melts at 
higher temperatures. 

We plot the (b(s)-b(O)) correlation function from the simulation in Fig. ||. For T > T x we obtain simple exponential 
decay but for T < T x we see an oscillation in the correlation function in agreement with Eq. ( p"3| ) . 

We plot the (t(s) • t(0)) correlation function from the simulation in Fig. We see the signature exponential decay 
of the correlation function of worm-like chains from which we can estimate the effective persistence length. The 
estimated persistence lengths are L p = 2.99 ± 0.01, 25.0 ± 0.005, 179.0 ± 0.001, correspondingly. 



VI. KINK-ROD STRUCTURE 

We show typical equilibrated conformations above, near, and below T x in Figs, pi ^|, |?]. The snapshots of the 
polymer configurations suggest that at low temperatures the polymer can be viewed as a collection of hard (straight) 
twisted rods that are connected by some kinks. This picture can be accounted for using a simple argument. We 
can model our system of two semiflexible polymers subject to the constraint of constant separation, as a semiflcxible 
ribbon, i.e. a semiflexible linear object with anisotropic rigidities whose Hamiltonian reads 



2 J ^ J V ds J ■ \ds 

where kij — «|| b^bj + kj_ (Sij — titj — bibj) determines the rigidity anisotropy of the ribbon, corresponding to bending 
parallel or perpendicular to the bond-director field. The ribbon structure would require Ku 3> (To be consistent 
with the hard constraint (see above) of constant separation of the polymers, we should take the limit of infinite /c||.) 
The partition function of a semiflexible ribbon in the — > oo limit can be written as 

2mb = >f m /m(.)Db(.)«p{-^ /d. (H) 2 + Ji, .b) 2 + *'[t,b]} (22) 

- / mwsbW^.bJexpj-^ / d, b]} , 

in which 7i'[t,b] controls the dynamics of b, and the functional delta-function enforces the constraint 

^ -1)00 = 0, (23) 
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to hold exactly at every point of the ribbon. Recalling dt/ds = H(s) n from the Frenet-Seret equations |18|, where 
H(s) is the curvature at each point and n is the unit normal vector to the curve, we can write the constraint as 



H(s) n(s) ■ b(s) = 0. (24) 

This constraint requires that at each point either H(s) = 0, which corresponds to a straight (rod-like) segment that 
can be twisted, or n(s) • b(s) = 0, which corresponds to a curved (kink-like) region where the the bond-director is 
locked in to the perpendicular direction to the curve normal, i.e., the binomial. 

We expect the (core) length of the kink regions to be very short at low temperatures, as observed in Fig. ^. We 
note that the conformational entropy of the chain is due to the degrees of freedom in the kink regions, whereas the 
twist entropy comes from the degrees of freedom in the rod segments. The average separation between neighbouring 
kinks is of the order of the persistence length. The ribbon thus tends to keep the rod segments as long as possible to 
maximally explore the twist degrees of freedom, while it can recover the same conformational entropy as a wormlikc 
chain from pivotal moves in the kink regions. This explains the kink-rod structure in low temperatures (a -C t p ). 
As the temperature increases, the kinks get closer to each other, until at some temperature their average separation 
becomes comparable to their size (£ p ~ a), and the kink-rod pattern disappears. 

This analysis can be understood in the context of the mean-field (b,t) model above. By observing that at low 
temperatures £bp <C £tp, one can imagine that there are roughly speaking rodlike (straight) segments of length £tp, 
each supporting a number of shorter segments of length £bp that are twisted, but de-correlated with one other. One 
can see that £bp is equal to the length scale A at which the strands undergo conformational fluctuations of the order 
of their separation: a 2 =< r 2 >= J 1 ^ x dq/£ p q 4 ~ X 3 /£ p . Hence, segments of length £bp are straight (for a <C £ p ). 
However, fluctuations of order a are sufficient to wash out the memory of twist. The anti-correlation in fact comes 
from the frustration as explained above. 

As the temperature is raised, the number of twisted rods in each segment N — £^-p / £bp decreases very quickly, 
until it saturates to unity at T — T x . For higher temperatures the mechanism changes, and the bond correlations 
are cut off by the tangent fluctuations. Hence, the short-range twist order does not survive anymore. All the main 
features of the above picture have also been observed in the simulation. 

A kink-rod structure similar to the one discussed here has indeed been observed in experiments done on actin 
filaments p9| . Actin is a charged polymer, and the mutual electrostatic repulsion of its different segments plays a 
major role in its structural stiffness. It is well known, however, that the introduction of multivalent counterions (ions 
of opposite charge that are necessary to neutralise the solution) can reduce the electrostatic repulsion, and even lead 
to attraction between like charged polymers, or different like charged segments of a same polymer ]2(J. In a recent 
experiment, Tang et al [Hj] used fluorescence microscopy techniques to image condensed (or collapsed) actin bundles 
that are formed due to the presence of multivalent counterions. Snapshots of the bundles, showing their typical 
conformations, are shown in Fig. [l^. A remarkable feature in the observation was the presence of sharp corners which 
connect relatively straight segments of the actin bundles, as can be seen in Fig. |l0|. Tang et al jl9| observed that this 
feature is present only when the bundle is made up of two or more filaments of actin, and is absent when there is only 
a single filament in the condensate. 

It is plausible to assume that the observed kink-rod structure can be accounted for by similar arguments to the 
one developed above. The only difference is the fact that the structure in the experimental case is a ring, as opposed 
to a chain with free ends. However, we do not expect this constraint to affect the argument, because the inherent 
competition between twist and bend degrees of freedom is local. Of course, more experimental efforts are needed 
to rule out other possible scenarios for the formation of the kinks, such as defects in the packing of more than one 
filament, sequence disorder, or metastable effects due to the dynamics of the collapse. 



VII. MEAN FIELD THEORY: CONFINED CHAINS 



In this section, we study the effect of confinement on double-stranded semiflexible polymers. We confine the double- 
stranded polymer in a d_L-dimensional subspace to a box of size R, while leaving it free in the remaining d\\ = d — d± 
dimensions. 

In terms of the tangent to the mid-curve t = dr/ds, and the bond-director b, the total Hamiltonian of the system 
can now be written as 



Hconf _ "H + "Hm 9_ 
k B T ~ k B T t% 



dsri, (25) 



where g is a dimensionless constant in addition to b, c, e defined in Section [II. The constants will be determined 
self-consistently by demanding that the relevant constraints hold on average: inextensibility of the individual strands, 
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constant separation between the chains, the bond-director being normal to the strands, and confinement of the polymer 
to a box of size R in d± dimensions ((r^_(s)) = R 2 ). This final constraint is valid as long as the chain length is much 
larger than R. We obtain the following set of equations for 6, c, e, and g: 



l<7. 



r-Rf , _d±_ _ 



(26) 



327? 4 



u = o 

Although the above nonlinear set of equations are very difficult to solve, we can get the behaviour of the solutions 
by looking at the asymptotics in the limiting cases, as summarised in Table |[ The full solutions are in fact smooth 
interpolations between the asymptotics. We solved the equations numerically for the experimentally relevant case; 
dx = 1 and d = 3. 

Having determined the constants self-consistently, we can calculate the correlation functions. For the tangent- 
tangent correlation in the parallel direction one obtains 



(t||( S )-ty(Q)> = 

whereas for the perpendicular direction one obtains 

d± 



W2b 



exp 



<t±(s)-t x (0)> = 



8^ 



2b 



exp 



(27) 



- (b- v / 6^ r g) 1/2 exp 
Similarly, for the bond-director field it yields 



1/2 



(*>(*)• b(0)) 



exp 



-{b-Vb^c) 1 ' 2 f- 



exp 



2a 2 Vh 2 



{b-Vtf^y/ 2 



(b + s/W^cfi 2 



(28) 



(29) 



while the rest of the two-point functions (the cross terms) are zero. The parallel component of the tangent-director 
correlation function decays purely exponentially. However, the correlation function of the perpendicular component 
of the tangent-director field, as well as that of the bond-director field, develop a crossover from purely exponential 
decay for b 2 > g and b 2 > c, to oscillatory decay for b 2 < g and b 2 < c, respectively. The phase diagram of the system 
in the space of dimensionless parameters R/£ p and a/£ p , is shown in Fig.^. The boundaries between different regions 
are obtained from solutions of Eq. ( p6| ) . 

It is instructive to examine the perpendicular component of the tangent-director correlation function in the limiting 
case b 2 <C g which corresponds to R/£ p <§C 1 (see Table 1). Using the asymptotic forms for b and g, one obtains 



(t_L( a )-t_L(0)> 



V8v2 



d 2 ± 



R 2 



1/3 



x exp 



\U P R 2 



1/3 



IU P R 2 



1/3 



(30) 



The effects of confinement are best seen in this limiting expression. The persistence length of the polymer in the 
confined directions is reduced to £± ~ (ipR 2 ) 1 ^ 3 . This 'deflection length' || is in fact the length at which roughening 
of a semiflexible chain of bare persistence length £ p becomes comparable to the confinement size (or separation of the 
confining walls), R: R 2 =< r 2 ^ >= fu i± dq/£ P q 4 - In other words, the presence of the boundaries provides another 
competing mechanism to cut off tangent correlations in the directions of confinement. Moreover, the oscillatory form 
of the correlation function with a period A = lj_ ~ (ipR 2 ) 1 ^ 3 implies a sinusoidal packing of the polymer in the 
confining cavity, where again the size of the oriented se gme nts (the period of the oscillations) is set by the walls 
cutting off the roughness of the polymer. As seen in Sec.[Vl[ a semiflexible ribbon develops a kink-rod structure at 
finite temperatures (due the strong anisotropy in rigidity), in which rod-like segments (about a persistence length 
long) are connected by rather sharp kinks with a core size of the order of the diameter of the ribbon. If such a 
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structure is restricted to a size less than the bare persistence length, new kinks have to be created at the confining 
walls to squeeze the ribbon into the available space. This therefore leads to the compact zigzag conformation of kinks 
and rods and the subsequent oscillatory tangent correlations. 

It is straight forward to calculate the free energy of the system. We find, 

Tt = ¥t + W V b + v 7 ^) 172 + Q>- V^) 1/2 + (b+ V^~9) 1/2 + (b- V^~g) l/2 \ ■ (31) 

Kb± KB-L Z£p L J 

It is interesting to examine the limiting behaviour of the free energy as a function of R. We obtain 

(d ± /2) 1/3 (e p /Rf /3 , for R/£ p < 1 



F(R) ^ L 
k B T ~ 2L, 



(a d ± /d 2 ) {Ip/Rf , for R/l v > 1 



(32) 



where a is a smoothly varying numerical coefficient ranging from a — 4 for a/£ p <C 1, to a = 1 for a/£ p 1. The 
free energy, interestingly, interpolates between the steric repulsion of non self-avoiding flexible polymers confined to 
R |21j, which has 1/R 2 behaviour, to Helfrich undulation free energy of stiff polymers confined to R H^2|, which has 
lfR 7 ? 3 behaviour (23). 

Using the above results, we now analyse some recent experiments by Ott et al p4J who measure the persistence 
length of actin filaments confined between microscope slides. The separation of the slides was of the order of R — Ifim 
and they found a persistence length (assuming a two dimensional worm-like chain) of L® = 16.7/im. These results are 
in the regime addressed by our model. For very small but finite R the chain fluctuates between the two plates. ^From 
the analogue of Eq. (|30|) calculated for single-stranded semiflexible polymers, we find (tj_) = (R/ Lp ^) 2 / 3 /2. Therefore 
one must include the fluctuations perpendicular to the confining plates in the calculation of the true persistence 
length. This implies that on average the polymer makes an angle 8 given by sin (9 = J {t 2 ) with the plates where 

9 < 1. One therefore has a corrected persistence length 4 rue « 4 0) /cos# = 4 0) /vi-<*i>- 

We therefore estimate 

for those experiments a correction of approximately 4%, i.e. Lp fue = 17.4/im. 



VIII. CONCLUSION 



In conclusion, we have calculated the properties of a well-defined model of a double-stranded semiflexible polymer 
and shown novel non-trivial differences between the high, low and zero temperature behaviour. At high T we find 
normal worm-like chain (WLC) behaviour and at low T we observe a novel kink-rod structure with short-range twist 
order whilst at T = we have a flat ribbon. 

In the analytical approach, the only approximation we have made is the relaxing of local constraints to global ones. 
Using a systematic l/d-expansion (see Appendix) we have shown that to calculate the 2-point correlation functions 
this is a valid approximation as higher order corrections only change the values of parameters but do not change the 
analytic form of the functions. Extensive MD/MC simulations confirm the analytical results. 

We have also examined the effect of confinement on the behaviour of semiflexible double-stranded polymers, and 
found four interesting regimes of the conformation and internal twist structure of these polymers, as summarised in 
Fig. O: (A) Weak confinement and relatively short bonds lead to free wormlike chain conformations with short-ranged 
twist anti-correlations, (B) weak confinement and relatively long bonds give rise to free wormlike chain conformations 
and twist disorder, (C) strong confinement and relatively short bonds yield sinusoidal packing of the chains and 
short-ranged twist anti-correlations, and finally (D) strong confinement and relatively long bonds lead to sinusoidal 
packing of the chains and twist disorder. 

There are a number of advantages evident in our approach. First, we introduce a microscopic model which remains 
true to the chemical structure of many biomolecules. Second, our approximate method of solving this model also lends 
itself to the analysis of the fluctuations in the system and to study intermediate-scale behaviour as well as the ground 
state (long length-scale) properties. Finally, this method could be easily extended to describe multi-stranded objects. 
We expect that the effect of an intrinsic twist changes the ground state but does not change any of the conclusions 
of our description, although we expect it to make the effective persistence length much higher. We hope to address 
such questions in a subsequent publication. 
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APPENDIX A: l/D-EXPANSION 



To justify the mean field approximations used in Sees. Ill and VII above, we perform a systematic 1/d-expansion 
that allows us to implement the constraints in a controlled way. Our approach is similar to the one successfully used 
by David and Guitter to study the crumpling transition of crystalline membranes |lj] . For simplicity, we consider a 
more primitive model of an elastic ribbon with the Hamiltonian 



ibbon 



A 



k B T 



dt 



- / As 



db 



-: t t • — 



d.s 



db 



As 



(Al) 



in which the tangent and the bond-director fields are subject to the following constraints:!; 2 = 1, b 2 = 1, t • b = 0, 
and b • (At/As) = (Eq.p3[)). A similar calculation with the more realistic (and more complicated) Hamiltonian of 
Eq.(|J) will lead to essentially the same conclusions. The partition function of the ribbon can be calculated as 



Z= VtVhVX 1 VX 2 VX z VX 4 e 



-S[t,b,Ai,A 2 ,A 3 ,A4] 



(A2) 



with 



S 



A 

+i 



As 
As 



C 
~2 



As 



Ah 

dl 



X 1 (s) (t 2 - 1) + A 2 (s) (b 2 - 1) + Aa(*) (t • b) + X 4 (s) (t • ^ 



(A3) 



in which {A a (s)} are the "stress" (Lagrange multiplier) fields enforcing the constraints. The integrations over t and 
b are now Gaussian, and can be performed to yield 



S c s[Xi] = - lndct 



-Ad 2 s + 2iXi iX 3 + i\ 4 d s 
iX 3 - id s X 4 -Cd 2 s + 2iX 2 



-id J As{Xi + A 2 ). 



(A4) 



Note that we have rescaled A, C, and {A a (s)} by d. Extremising the effective action corresponds to the saddle point 
or mean field solution (d = oo). The saddle point equations yield iXi = 1/(8A), iX 2 = /(8C), and A3 = A4 = 0. 
To proceed to the higher orders in 1 /d, we need to calculate the A-propagators defined as 



(X a (s)Xp(s')) = 



s 2 s cS 



A straightforward calculation then leads to 



Ai(«)Ai(- 


-q) 


X2(q)~X 2 (- 


-i) 


X 3 (q)X 3 (- 


-Q) 


X 4 (q)X 4 (- 


-q) 


A 3 (<Z)A 4 (- 


-q) 



5X a (s)5Xp(s') 



saddle— point 



A 2 
= 8d q 

C 2 



A 



0(1), 
0(1), 
0(1), 

0(l/q 2 ), 



2d 

A + C 2 



(A5) 

(A6) 
(A7) 
(A8) 
(A9) 
(A10) 



while all the others are zero. Note that we have only kept the large momentum limit, since we are interested in the 
local (short distance) behaviour of the Lagrange multipliers |l4| . 

The above A-propagators, the correlators for the tangent and the bond fields 
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and the three-point vertices 



(ii(qM-q)) 

Hq)bj(-q) 



f + \/(AAy 

Si 



'1.1 



Cq 2 + l/{ACy 



Uiqfoiq'^i-q-q') 

Uq)bj(q')H-q-q') 

ti{q)bj(q')Xs(-q-q') 
Uiqfbjiq'yXii-q - q 1 ) 



2i5i 



2i5i 



= i5, 



-q'Si 



(All) 
(A12) 

(A13) 

(A14) 

(A15) 

(A16) 
(A17) 



as read from Eq.(A3), could now be used to construct diagrammatic expansions. Examining the 2-point correlation 
functions for t and b, we then find that the perturbative expansions are well-behaved (not singular), and only correct 
the numerical values of the coupling constants by finite amounts at each order. We thus conclude that the mean field 
behaviour corresponding to the saddle point approximation is qualitatively valid. 
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FIG. 1. The schematic of the double-stranded semiflexible polymer of two chains separated by a distance a. Note the bond 
director field b(s). 




FIG. 2. The solution of the self consistent equations for the constants b and c as a function of u — a/£ p in d — 3. The value 
u c ~ 4.27/d corresponds to the transition point. 
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FIG. 3. Coarse-grained model - Like and Unlike twists meeting. The bond vector B is rotating about the tangent vector in 
the same direction for like twists and opposite direction for unlike twists. 
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FIG. 4. The schematic of the double-strand model used in the simulation. 
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FIG. 5. Typical conformations from MD/MC simulations of a ribbon made up of two chains of 400 monomers above T x 
k h = 1. 
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FIG. 6. Typical conformations from MD/MC simulations of a ribbon made up of two chains of 400 monomers near T : 
k b = 10. 
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FIG. 8. The (b(s) • b(0)) correlation function measured in the simulations for temperatures kt = 1,5,50 corresponding to 
b 2 > c,b 2 ~ c and b 2 < c. The averages were done over ~ 10 4 statistically independent samples. The error bars are the size of 
the symbols. 
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s 

FIG. 9. The (t(s) ■ t(0)) correlation function measured in the simulations for temperatures kt — 1, 10, 50 corresponding to 
b 2 > c, b 2 ~ c and b 2 < c. The averages were done over ~ 10 4 statistically independent samples. The error bars are the size of 
the symbols. 

FIG. 10. Snapshots of condensed actin rings due to the presence of multivalent counterions, from the experiment by Tang 
et al The kink- rod structure is manifest, in bundles with two or more filaments [H. 
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TABLE I. The limiting behaviour of the confined polymer: Asymptotic behaviour of the solutions of Eq.(B 
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Structural Variants of F-actin Bundles 




